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ABSTRACT 

We present the results of SCUBA-2 observations at 450 fim and 850 fim of the field lensed by the 
massive cluster A370. With a total survey area > 100 arcmin^ and 1 a sensitivities of 3.92 and 0.82 
mJy/beam at 450 and 850 iim respectively, we find a secure sample of 20 sources at 450 fj,m and 26 
sources at 850 /im with a signal-to-noise ratio > 4. Using the latest lensing model of A370 and Monte 
Carlo simulations, we derive the number counts at both wavelengths. The 450 fim number counts 
probe a factor of four deeper than the counts recently obtained from the Herschel Space Telescope 
at similar wavelengths, and we estimate that ^47-61% of the 450 /im extragalactic background light 
(EBL) resolved into individual sources with 450 /im fluxes greater than 4.5 mJy. The faint 450 /im 
sources in the 4 a sample have positional accuracies of 3 arcseconds, while brighter sources (signal- 
to-noise > 6a) are good to 1.4 arcseconds. Using a deep radio map {la ^ 6 /iJy) we find that the 
percentage of submillimeter sources having secure radio counterparts is 85% for 450 fim sources with 
intrinsic fluxes > 6mJy and 67% for 850 /im sources with intrinsic fluxes > 4mJy. We also find that 
67% of the > Aa 450 ^m sources are detected at 850 /im, while the recovery rate at 450 /im of > 4 cr 
850 /^m sources is 54%. Combined with the source redshifts estimated using millimetric fiux ratios, 
the recovered rate is consistent with the scenario where both 450 //m and 20 cm emission preferentially 
select lower redshift dusty sources, while 850 /zm emission traces a higher fraction of dusty sources at 
higher redshifts. We identify potential counterparts in various wavelengths from X-ray to mid-infrared 
and measure the multiwavelength photometry, which we then use to analyze the characteristics of the 
sources. We find three X-ray counterparts to our robust submillimeter sample (S/N > 5), giving an 
active galactic nucleus fraction for our 450 (850) /im sample of 3/8 (3/9) or 38% (33%). We also find 
a correlation between the Kg-hand and the 850/im/20 cm flux ratio. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: formation — submil- 
limeter 
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1. INTRODUCTION 

Based on studies of the integrated light originating 
from outside the Milky Way galaxy — the extragalactic 
background light (EBL) — we now know that there is a 
comparable amount of light absorbed by dust and re- 
radiated in the far-infrared (FIR) /submillimeter (100 fim 
< X < 1000 /im) as there is seen directly in the opti- 
cal/U V ( |Puget et al.|1996[|Fixsen et aL|1998"||Dole et al. 
20061, which indicates that much ot the star formation 
and active galactic nucleus (AGN) activity in the Uni- 
verse is hidden by dust. Understanding dusty star forma- 
tion is therefore critical to a full understanding of galaxy 
formation. 

The first step in such work is identifying the individ- 
ual dusty galaxies. Observations have been made to- 
wards this goal using both the small space-ba sed Her- 
schel Space Observatory (hereafter Herschel; Pilbratt 



et al.||2010[ ) and, at longer wavelengths, ground- based 
telescopes (e.g., the James Clerk Maxwell Telescope 
(JCMT)). However, the fraction of sources that can be 
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identified with such observations has a fundamental limit 
due to the poo r resolution. This limit is known as th e 
confusion limit ( |Scheuer|195"7||Condon|f974{ |Hogg|200T| ) , 
which is the maximum sensitivity an observation is able 
to reach for a given beam size due to the overlap of 
sources on the map. Recent observations carried out by 
Herschel at 250, 350, and 500 /im can only resolve a small 



amount of the EBL (15, 10, and 6%; Oliver et al. 

With a 15 meter dis' 
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due to the confusion limit 
JCMT blank-field su bmillimeter surveys with SCUBA 
([Holland et al.|[l999| ) only resolved ~20%-30% of the 
850 /im EBL into distinct, bright submillimeter galaxies 
(SMGs) with Ssso^m > 2mJy (e.g., [B arger et al. 1998 

" - - lyyg 



Hughes et al.. 1998, Barger et al. 



|1999aj jEales et al. 



'Coppin et al. 2006) before reaching the co nfusion limit 
'i'he most recent surveys using LABOCA ( [Siringo et al. 
2009[ ) on the APEX telescope at 870 /im obtained similar 
results (Weifi et al. 2009). 

The most fundamental consequence of the confusion 
limit is that the poor resolution at 850 /im prevents the 
study of faint SMGs below 2 mJy in blank fields. A 
small number of fainter SMGs with intrinsic fiuxes be- 
tween 0.1 and 2 mJy have been detected from observa- 
tions of the fields of massive lensing clusters. In these 
fields, the intrinsically faint fluxes of background sources 
are gravitationally amplified to a detectable level. The 
sources found in t his way contribute more than 50% of 



the 850 /im EBL (Small et al. 1997 2002 Cowie et al 



2002 Knudsen et al.||20D5 |. However, the positional un- 



2 



certainties of the 850 /zm sources can cause large uncer- 
tainties in the amphfications and in the intrinsic source 



fluxes (Chen et al. 2011). It is therefore essential to con- 
duct surveys on lensmg cluster fields with as high a spa- 
tial resolution as possible to circumvent the confusion 
limit and to construct the most accurate possible num- 
ber counts distribution. 

The poor resolution (e.g., 14" FWHM on the JCMT at 
850 /im) also makes identifying the counterparts in other 
wavelengths difficult and time consuming. Deep, high- 
resolution radio, mid-infrared (MIR), optical, UV, and 
X-ray observations have been used to identify candidate 



tion and active galactic nucleus (AGN) activity (Barger 


et al.| 2000! 'Chapman ct al."2005| Pope 


et al. 200B Ivi- 


son et al.|f2007; Gcorgantopoulos et al.| 


201 ip. However, 



servations have shown that due to their clustered nature 



( [Weifi et al.||2009[ ), SMGs are likely to be resolved into 
multiple sources and that different tracers identify differ- 
ent counterparts based on their own biases (e.g. , Wang 



et al.l2011| [Barger et al.pOT2l [SSIoTcic et al.|2012| [K arim 



et al.||2012[ ). 

Many of these problems can be avoided by observing 
at 450 /xm rather than at 850 fim. The higher spatial 
resolution at 450 /xm means that we can detect several 
times more sources before hitting the confusion limit, and 
the more accurate positions allow us to better determine 
the amplifications when observing lensed sources. In ad- 
dition, while 850 /im observations have the unique ad- 
vantage of tracing the star-forming galaxies to very high 
redshifts (z ~ 10), thanks to the negative if-correction 



(Blain et al. 2002), it has been argued that owing to 
the location on the Rayleigh- Jeans tail of the blackbody 
spectral energy distribution (SED), the 850 /xm sel ection 
is biased against high dust temperature sources ( Blain 



et al. 2004 Chapman et al. 



Casey et al. 2005) 



'lb mitigate this selection ettect, observations at rest 
frame wavelengths close to the peak of the blackbody 
SED (~ 100 - 200 ^m) that are less affected by the 
dust temperature are needed. The SPIRE instrument 
mounted on Herschel is designed to do this job; however, 
its ability to detect sources is again heavily restricted by 
the confusion limit due to the small aperture size of the 
telescope. 

450 UTa mapping was very inefficient with SCUBA 
( Chapman et al. 2002 Small et al. 2002 1. Based on a 
handful of extracted sources (< 10) with 450 fiui fluxes 
above 10 mJy, only 15% of the 450 fim EBL was resolved 
dSmail et al!][2002| . The advent of SCUBA-2 ( [Holland 
et al. 2006| on the JCMT with its two orders of mag- 
nitude more pixels/bolometers than on SCUBA greatly 
boosts our ability to resolve the bulk of the 450 /im EBL 
and to detect galaxies not biased by dust temperature. 
SCUBA-2 observes simultaneously at 450 and 850 /im via 
a dichroic beamsplitter. Unlike SCUBA, the bolome- 
ters on SCUBA-2 are optimized for observing at 450 and 
850 /im, and the filter transmission, especially at 450 /im, 
is much better than that of SCUBA. Moreover, the field- 
of-view of SCUBA-2 is more than 10 times larger than 
that of SCUBA. With all of these new features, SCUBA- 
2 is an order of magnitude faster than SCUBA, which 
makes it possible to finally study the 450 population 
in detail. The amount of FIR EBL that can be resolved 



by SCUBA-2 is unprecedented as we show later in the pa- 
per, and demonstrates both the power of SCUBA-2 and 
also the promise of the next generation of submillimeter 
telescopes (CCAT). 

Here we describe the results that we have obtained on 
the lensing cluster field Abell 370 (A370) from our on- 
going SCUBA-2 program to survey a sample of massive 
galaxy clusters at 450 /im. The goals of this program are 
to resolve the 450 /im EBL as fully as possible, to effi- 
ciently identify the counterparts of the 850 /im sources, 
and, most important of all, to study the nature of the 
SMC population in as unbiased a fashion as possible. In 
Section 2, we discuss our SCUBA-2 observations, data 
reduction, and ancillary data. Since the results involve 
a relatively new instrument we describe the methodol- 
ogy in detail. Readers who are primarily interested in 
the results can skip directly to section 3. In Section 3, 
we describe our source extraction and present our cata- 
log. In Section 4, we determine the number counts and 
the amount of the EBL that our data resolve at both 450 
and 850 /im. In Section 5, we discuss the properties of the 
sources at submillimeter and radio wavelengths. Using 
the rich set of ancillary data on this field from the radio 
to the X-ray, in Section 6, we analyze the properties of 
the counterparts to the submillimeter sources, construct 
color-color and color-flux diagrams, and examine the na- 
ture of both 450 and 850 /im sources. Finally, we summa- 
rize our results in Section 7. Throughout the paper, we 
adopt the AB magnitude system for the photometry. We 
also assume the Wilkinson Microwave Anisotropy Probe 
cosmology of Hp = 70.5 km s~ ^ Mpc^^, flM = 0.27, and 
r^A = 0.73 ( [Larson et aLl|2011[ ). 



2. OBSERVATIONS, DATA REDUCTION AND ANCILLARY 

DATA 

We carried out the SCUBA-2 observations in Decem- 
ber 2011 and in February and March 2012 (Program IDs: 
M11BH26A and M12AH26A) under the best weather 
conditions (band 1, tau225 GHz < 0.05). A370 has a large 
Einstein ring with an angular diameter size of 3'. To ex- 
ploit the strong magnification and to efficiently use any 
available time to discover faint sources, we adopted the 
CV Daisy patterrj^ to perform the observations. This 
observing mode is optimal for an area less than 4'. Each 
scan/cycle takes around 20 minutes to complete. Flat- 
fields were taken at the beginning and end of each scan, 
however only the first one is used for data reduction. 
The extinction correction coefficients are obtained from 
the JCMT water-vapour radiometer (WVM) using the 
183 GHz water line to calculate the precipitable water 
vapour (PWV) along the line-of-sight of the telescope 
( [Dempsey et aL[[2012p . 

We reduced the data using the Dynamic Iterative Map- 
Maker (DIMM) in the SMURF package from the STAR- 
LIN K software developed by the Joint Astronomy Cen- 
tre ([Jenness et al. 2011) released before July 1, 2012. 
The program automatically performs dark subtraction, 
flatfielding, extinction corrections, and data fiagging. It 
also removes the DC steps from each bolometer, performs 
frequency domain filtering, and concatenates the reduced 

^ A detailed descriptio n of the CV Daisy pattern can be 
found at http : / /www . j ach .hawa il ■ edu/ JCMT/continuuin/scuba2/| 
scuba2_obsmodes .html#scan 



Fig. 1. — Demonstration 850 /im maps for post-processing, (a) Map created from DIMM, (b) Map in (a) convolved with a broad Gaussian 
(FWHM = 30" for 850 fim). (c) Map in (a) subtracted from map in (b). (d) Result of performing matched- filter on map in (c). Note that 
these are all S/N maps scaled to maximize the visibility of the sources. 



data from each chip to produce the final sky map. DIMM 
provides considerable freedom to the user to fine-tune the 
parameters that control the jobs mentioned above, and 
those parameters are stored in a configuration file that 
is included in DIMM during the data reduction. 

There are four detector chips in each waveband with 
1280 (40x32) bolometers per chip. We ran DIMM on 
the individual chips in each scan and used the MO- 
SAIC JCMTJMAGES recipe in PICAR D (Pipeline for 
Comb ining and Analyzing Reduced Data; IJenness et al.| 
2008| to mosaic and coadd the maps. We found this 



method produces better signal-to-noise (S/N) maps than 
running DIMM on all the data at once. The noise maps 
and the exposure time maps were also created from 
DIMM. The noise map is obtained by computing the 
variance of the data that lands in each pixel. We tested 
the robustness of the noise map created from DIMM by 
checking the standard deviation of the S/N maps. If the 
noise map is robust, then the standard deviation of the 
S/N map should be close to 1. All the noise maps are 
accurate with an underestimation less than 5%. We sus- 
pect that the cause of the small underestimation is the 
correlated noise from the large scale structure caused by 
the atmospheric noise (see Figure [TJa)), which is ulti- 
mately subtracted out in our post-processing technique 
(Section 2.1). 

2.1. Flux Calibration 

We adopted the configuration file dimmcon- 
fig_bright_coTnpact.lis provided by SMURF for DIMM 
to produce the calibrators' maps in units of picowatts 
(pW), and we used the CALC_SCUBA2_FCF recipe 
from PICARD to calculate the Flux Conversion Factors 
(FCFs) for the SCUBA-2 maps. 

The recipe performs a Gaussian fit for the calibrator 
and computes the total flux in units of pW arcsecond^. 
It then divides the expected flux in Jy by the com- 
puted flux to obtain the FCF_[ARCSEC] in units 
of Jy/pW/arcsecond^. It has been reported that 
FCF_[ARCSEC] is very stable from extensive exam- 
ination of the calibrators observed over the period 
of SCUBA-2 commi ssioning and science verification 
( Dempsey et al.|[2012 ). We doubled checked this on our 
own calibrators and confirmed that the number is indeed 
stable. 

The majority of SMGs are found to be located at high 
redshifts {z > 1), with the median z ~ 2.5 ( |Chap-, 



|man et al.||2005[ [Wardlow et al.||2011|. Based on high- 
resolution interferomctric observations, a physical reso- 
lution l ess than a few kpc is needed to resolve individual 
SMGs dKnudsen et al.|2010[ ). With our 7'.'5 resolution at 
450 /xm, we can barely resolve sources at z '--^ 0.1. Thus, 
it is safe to assume that SMGs are mostly unresolved 
point sources at the resolution of our maps. 

To optimize the point source extraction, we adopted 
the matched-filter algorithm, which is a maximum 
likefihood estimator of the source strength, for post- 
processing in o rder to increase the S/N for point source 
detection (e.g., 'Serjeant et al. 2003). Assuming S(i,j) and 
cr(i,j) are the signal and r.m.s noise maps produced by 
DIMM and PSF(i,j) is the signal point spread function, 
the filtered signal map F(i,j) would be 



J:^,^S{^,J)/ai^,JrxPSF{^,J)] 



E,,,[l/a(z,jTxP5F(z,j7] 



and the filtered noise map N(i,j) would be 

_/V(i j) = - 



(1) 



(2) 



Ideally, the PSF for the matched-filter algorithm is a 
Gaussian normalized to a peak of unity with FWHM 
equal to the JCMT beamsize at a given wavelength (i.e., 
7'.'5 at 450 /im and 14" at 850 /im). However, the map 
produced from DIMM usually has low spatial frequency 
structures that need to be subtracted off before perform- 
ing the source extraction. Thus, before running matched- 
filter, we convolved the map with a broad Gaussian (22" 
at both 450 and 850 /Ltm) normalized to a sum of unity 
and subtracted this convolved map from the original 
map. Note that we experimented with different FWHM 
sizes of the broad Gaussian used to convolve the maps for 
background subtraction, and we found that the source 
fluxes and the S/N are not sensitive to the size of the 
FWHM (with variations lower than 5%) for reasonable 
choices. We adopted 22", which gives a good suppres- 
sion of the negative signals. In Figure [l] we show the be- 
fore and after versions of the background subtraction and 
matched-fllter. To determine the fluxes, we processed 
the PSF used for matched-filter similarly. It becomes a 
Gaussian with a convolved broader Gaussian subtracted 
off, which gives a Mexican hat-like wavelet. We adopted 
the PICARD recipe SCUBA2_MATCHED JILTER for 
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TABLE 1 

FWHM OF THE PSFS AND FCFS OF THE CALIBRATORS 



Date 


Calibrators 


FWHMx 


FWHMy 


FCFJBEAMMATCH] 






(arcsec; 450/850) 


(arcsec; 450/850) 


(Jy/Bcam/pW; 450/850) 


122611 


CRL618 


8.96 / 13.26 


8.64 / 13.03 


831.513 / 606.928 


122811 


Uranus 


8.74 / 13.06 


8.41 / 13.23 


750.590 / 567.313 




Uranus 


9.41 / 13.40 


9.51 / 13.94 


819.377 / 602.322 


020412 


Uranus 


8.32 / 12.77 


9.44 / 13.65 


697.301 / 568.398 


020512 


Uranus 


8.47 / 13.06 


8.71 / 13.45 


710.364 / 563.341 


Average 




8.78 / 13.11 


8.94 / 13.46 


761.829 / 581.660 


Stddev 




0.43 / 0.24 


0.50 / 0.36 


61.454 / 21.111 



TABLE 2 

SCUBA-2 OBSERVATIONS ON A370 FIELD 

R.A. (J2000; h:m:s) 02:39:53 

DEC. (J2000; °:':") -01:34:38 

Area (850 /im; arcmin^)" 124.579 

Area (450 /im; arcmin^)" 121.095 

Central Noise (850 /im; 1 tr; mjy/beam) 0.82 

Central Noise (450 ^im; 1 a; mJy/beam) 3.92 

" Total area to 3 times the central noise level 



the tasks described above. 

Because of the atmospheric fluctuation, especially at 
450 /Ltm, the signal distribution is not always a perfect 
Gaussian with the expected FWHM, which could lead 
to incorrect flux measurements if we adopt the stan- 
dard PSF for the matched-filter. However, this un- 
certainty can be taken into account by going through 
the same source extraction process on the calibrators 
as was done on the science maps. We therefore applied 
matched-filter to the calibrators with the same PSF func- 
tion used on the science maps, including the background 
subtraction, and divided the calibrators' expected flux 
by the after-matched-filter fitted peak flux to obtain 
FCF_[BEAMMATCH] in units of Jy/beam/pW. This is 
the FCF that we used to calibrate our science maps. 

There are several calibrators available each night, and 
the variation of FCF_[BEAMMATCH] between cahbra- 
tors is usually larger than the statistical error. We chose 
the primary calibrator Uranus most of the time, because 
of its high S/N. Moreover, compared to the other primary 
calibrator Mars, the angular size of Uranus is more sta- 
ble and smaller than the beam size (unresolved), which 
is ideal for flux calibration for point sources. However, 
in cases where the phase of Uranus is highly distorted, 
we considered the secondary calibrator. Table [T] shows 
the 2-D gaussian fitted FWHM on the x and y axes, 
along with the FCF_[BEAMMATCH] of aU our cahbra- 
tors. The FWHM of the PSFs and the FCFs are very 
stable for all our observations (less than 10% variations), 
which is evidence that the PSFs at both wavelengths are 
well determined. 

2.2. Science Maps 

We adopted the configuration file dimmcon- 
fig_blank_field.lis for our science maps. The biggest 
difference between dimmconfig-blank_field.lis and dimm- 
config -bright-Compact. lis used for calibrators is that the 
latter is less aggressive in identifying bad bolometers 
and DC steps, in order to ease the misidentification 



between bad bolometers/DC steps and bright sources. 

For each night of data, we used DIMM to reduce the 
data on the individual chips in each scan and then coad- 
ded them to produce one science map and one calibra- 
tor map. We performed the matched-filter algorithm on 
both maps and applied the FCF_[BEAMMATCH] ob- 
tained from the calibrator to the science map. Finally, 
we coadded all the science maps that were created on a 
night-by-night basis to produce the final map for source 
extraction. 

Because we subtracted the original signal map from 
a self-convolved map during the post-processing to get 
rid of the large scale structure, negative troughs are gen- 
erated around strong sources (see Figure 1(d)). Those 
negative troughs may suppress nearby weak sources. 
We therefore re-did the post-processing by masking out 
sources with S /N greater than 4 and iterated this proce- 
dure until there were no new sources with S/N greater 
than 4. We then put the masked sources back into the 
masked images to generate a hybrid map, which is the 
final map that we used for our analysis. 

The noise map created through Equation 2 would in- 
evitably underestimate the noise level, because of the 
spatially correlated noise (Figure [l|[d)). We corrected 
for this underestimation by multiplying the noise map by 
the standard deviation of the original S/N map (Equa- 
tion 1 divided by Equation 2). In Figure |2j we plot the 
histograms of the final data maps at botn wavelengths 
within the effective regions. The effective regions are 
the regions where the noise level is less than three times 
the central r.m.s., meaning the region where the effective 
exposure time is more than 1/9 of the central exposure 
time. The black curves represent the region of noise. The 
excess negative signals at 850 /im come from one nega- 
tive source (see discussion in Section 3.3). We present 
the final central noise level, along with other observa- 
tional information on A370, in Table [2j 

2.3. Astrometry Calibration 
SMGs are often traced by 20 cm sources, thanks to 



et al.| 


1985 


Condon 


1992 


must 


3reak down at very hi 



Although this correlation 
redshifts due to the Comp- 
ton cooling between the relativistic electrons and the cos- 
mic microwave background photons, so far no clear ev- 
idence for_eTOllitiono£_this correlation has been seen to 
2 ~ 5 ( [Barger et aL]|2012 |. 

Since most SMGs found so far are located at z < 5, 
the radio sources detected in high-resolution maps taken 
by interferometers are a good sample to use for our abso- 
lute astrometry calibration. We searched for the offset of 
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the submillimeter maps that maximized the stacke d sub- 
millimeter s iRiial over all the radio sources in the Wold 
et al.| (20121 radio catalog of the field. We found a sys- 
tematicsHift of (0", 2") in our 850 /xm map and (-1", 1") 
in our 450 /im map. We corrected this small offset before 
proceeding to the source extraction. 

2.4. Ancillary Data 

A370 is one of the most extensively studied gravita- 
tional lensing cluster fields, with its multiwavelength ob- 
servations from the X-ray to the radio. In Table Isj we 
list all the data on A370 that we used in this paper, along 
with the corresponding references. 

The Chandra X-ray survey of the A370 field with 
66.6 ks of useful exposure time has a la sensitivity 
of ^ lO"""^^ erg cm^ s~^ in the 2-7 keV ban d and 
~ 3x 10"^^ erg cm^ s~\in th e 0.5-2 keV band ( |Bautz 
al][2000{ [Barger et al.l|2001[ ). 

i'he optical z-band data were taken by |Hu et a~ 
( 2010 1 using Suprime-Cam on Subaru for a study of high- 



A370 Optical, NIR, 



TABLE 3 
MIR, FIR, 



Radio Photometry 



redshift Lya emitters. The 5 a sensitivity for a 3" diam- 
eter aperture is 25.4 AB magnitude. 

The deep NIR J, H, Kg observations were carried out 
by Keenan et al. (2010) using ULBCam on the University 
of Hawaii 2.2 m (im2.2m) and WFCam on UKIRT. The 
5tT detection limits are 23.5, 23.0, and 23.0 AB magni- 
tude, respectively. We adopted their configu ration files 
and used SExtractor ( Bertin fc Arnouts|1996 ) to extract 
sources. 



The A37 field was first surveyed by Spitzer (Werner 
et aL|2004| ) in 2006 (PI: Giova nni Fazio; PIP: 137 ) in all 
the ik Arr, 



Xrray Camera (IRAC; Fazio et al. 2004) bands 
(3.6-8.0 /xm). Data with wider coverage at 3.6 and 4.5 /im 
were later taken in 2009 and 2010 (PI: Eiichi Egami; 
PID: 60034) when Spitzer entered the warm phase. We 
again used SExtractor to extract sources from the Spitzer 
maps. We measured the source fluxes using apertures 
with diameters of 4'.'8 for 3.6 //m and 4.5 /zm and 6'.'0 for 
5.8 and 8.0 /zm, which are roughly three times the 
FWHM of the PSFs and are the same as those used in 



IWang et al. (20041). We estimated the sensitivities using 
^ ' ^^The 



Gaussian fits to the fluxes measured at random source- 
free positions, finding 5cr limits of 23.3 at 3.6 /im and 
4.5 /im, 23.0 at 5.8 /fm, and 22.6 at 8.0 //m, aU in AB 
magnitudes. 

The deep VLA radio data of the field were taken using 
both the A and B configurations. The final image reaches 
a 5 (T limit of 28.5 /iJy /beam near the ce ntral region with 
a beam size of ~ 1'.'7 ( |Wold et al.|[2012 ). 

3. SOURCE IDENTIFICATION AND CATALOG 

3.1. Source Extraction 

On both maps, we first extracted sources that have a 
peak S/N greater than 3.0 within the effective regions. 
The way we extracted the source is that we found the 
maximum pixel within the desired region, took the po- 
sition and the value information of the peak, and sub- 
tracted a PSF centered at the source and scaled accord- 
ingly. We iterated this process until the peak S/N value 
was less than 3.0. We treated these peaks as the pre- 
liminary catalogs at 450 /im and 850 /xm. We then cross- 
correlated the two catalogs to find the counterparts in 
the other bandpass. We considered a counterpart re- 
covered if the position of the 450 /im source lay within 



Wavebands 



Instruments /Telescopes References 



X-ray (2-7 keV) 

X-ray (0.5-2 keV) 

z 

J 

H 

Ks 

3.6 

4.5 

5.8 

8.0 ^tm 
450/850 Atm 
20 cm 



ACIS-S3/ Chandra 
ACIS-S3/ Chandra 
Suprime-Cam/Subaru 
ULBCam/UH2.2m 
ULBCam/UH2.2m 
WFCAM/UKIRT 
IRAC/ Spitzer 
IRAC/ Spitzer 
IRAC/ Spitzer 
IRAC/ Spitzer 
SCUBA-2/JCMT 
VLA 



2010 
2010 



Barger et al. 2001 
Barger et al. 2001 
Hu et al. 2010 
Kccnan et al. 2010 
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Fig. 2. — The histograms of the S/N values of the pixels located 
within the regions where the noise level is less than three times the 
central r.m.s. (gray area). The black curves represent the region 
of noise. The excess negative signals at 850 /xm come from one 
negative source (see discussion in Section 3.3). 



the 850 /im beam. The only exception was A370-450.9, 
which lies slightly outside the beam area of A370-850.5 
but likely contributes a certain amount of the 850 /im flux 
and makes A370-850.5 appear elongated. We considered 
both A370-450.9 and A370-450.il to be counterparts to 
A370-850.5. 

To form the final catalogs, we kept every > 4 cr source 
in the preliminary catalogs, as well as every > 3 cr source 
that had a > 3 cr counterpart in the other bandpass. In 
Table [I] ([s]), we give the > 4 cr 450 (850) /im sources. 



6 



as well as the > 3 cr 450 (850) sources with > 3 a 
850 (450) /^m counterparts. In Column 1, we give the 
name of the source; in Column 2, our identification for 
the source; in Columns 3 and 4, the right ascension and 
declination in J2000 coordinates; in Column 5, the S/N 
in the discovery bandpass; in Column 6, the flux in the 
discovery bandpass; in Column 7, our identification for 
the counterpart in the other bandpass, where available; 
in Column 8, the S/N for that counterpart or the S/N 
measured at the peak position in the other bandpass; 
in Column 9, the flux for that counterpart or the flux 
measured at the peak position in the other bandpass. 

We give our assessment of the reliability of our source 
extraction in Section 3.2. We show our final S/N im- 
ages in Figure |3] overplotted with all the sources in Ta- 
bles |4] andjs] (red circles and blue squares for 850 /im and 
450 /im, respectively) and with the noise contours. In to- 
tal, we detected 26 sources at 850 /im and 20 sources at 
450 /fm. The number density is 3 (10) sources per 1000 
beam areas at 450 (850) /im, which indicates that our 
maps are not yet reaching the confusion limit (roughly 1 
source per 30 beam areas). Thus, the flux boosting and 
the positional errors caus ed by source confusion should 

Note that the maps are all 
We experimented with 



be negligible (Hogg 2001) 



gridded into 1" by 1" pixels, 
different pixel sizes and found that the results are not 
affected by the different choices. This is because the 
matched-filter technique is not affected by the pixel size 
(Equations 1 and 2), as long as the PSF is well sampled. 

3.2. Reliability of Source Extraction 

We used two methods to quantify the number of spu- 
rious sources above 4 a that might be included in our 
catalogs. First, we created a map with the same spa- 
tial dimensions as the real maps, populated the pixels 
with random Gaussian signal convolved with the PSF at 
the given wavelength, and then ran the peak identifica- 
tion. We iterated this process 100 times at each wave- 
length and found that we should expect 1.8±1.3 sources 
at 450 /im and 0.6±0.8 sources at 850 /im to be spurious. 

Second, we performed our same source extraction 
method on the inverted maps, where the negative sig- 
nal becomes positive and vice versa, to see how many 
negative sources are present at > 4 cr. We found 4 and 1 
at 450 and 850 /im, respectively, which agrees reasonably 
well with our first approach and means that we could 
have that many spurious sources in our > 4 cr catalogs. 

For the 3 cr sources identified as being present in both 
bandpasses, we again created random noise maps as de- 
scribed above, extracted sources between 3 and 4 cr, and 
then cross-correlated them with the > 3 cr sources ex- 
tracted from the real maps at the other bandpass. We 
iterated this process 1000 times and found 0.2 and 0.3 
spurious sources at 450 and 850 /im. In other words, 
sources identified at the 3 cr level at both wavelengths 
are even more reliable than > 4 cr sources only detected 
in one bandpass. 

3.3. Comparison with Previous Catalogs 

A370 was surveyed by SCUBA at both 450 and 850 //m 
within the circular region c entered at the cluster centroid 
with a radius less than 1^5 ( Cowie et al.|2002 Small et al.' 
2002| ). A370-850.3 was also observed ( |Chen et al..2011| ) 



with the hig h-resolution Submillimeter Array (SMA; Ho 
et all 12004 1. We list our detected sources within the 
8CUBA observed region along with previous results in 
Table [6) Our maps do not recover all the sources de- 
tected in |Cowie et al. (2002), because their maps have 
better sensitivity (< (J. 8 mJy/beam) at 850 /im. How- 
ever, most of the recovered sources agree with the litera- 
ture in terms of their fluxes and positions, exc e pt for the 



flux on A370 -450.1 given by Smail et al. (2002). Knudsen 



et al. (2008) found that it is diflicult to reliably calibrate 



the yCUBA 450 /im sources, because the PSF is not well 
described by a 2D Gaussian and is very sensitive to the 
deformation of the JCMT dish. Thus, the 450 /im fluxes 
obtained from SCUBA could be problematic. This prob- 
lem is significantly less in the SCUBA-2 maps, since, as 
we showed in Table IT] the PSFs at both 450 /xm and 
850 /im are stable andwell modeled. 

3.4. Completeness Tests 

We took the noise maps created in the previous sections 
and ran a completeness test. We randomly populated 
the image with sources of a given flux and then ran the 
source extraction on the maps down to 4 cr to see if we 
could recover the added sources. A source is considered 
recovered if it is detected above 4 cr and if its position 
is within the beam area. We iterated this process 5000 
times on each map with a flux step of 1 mJy (0.5 mJy) 
from 0.1 to 30.1 mJy for 450 (850) /im to estimate the 
completeness. We plot the results in Figure [4| The 50% 
completeness is around 4 and 17 mJy at 850 and 450 /im, 
respectively, and the 80% completeness is around 5 and 
21 mJy, respectively. 

4. NUMBER COUNTS 

A370 is one of the most extensively studied gravita- 
tional lensing cluster fields with a well-constrained lens- 
ing model. We used LENSTOOL, a software package 
that models the effects of gravitational lensing by taking 
three-dimensi onal mass distribut ions within the cluster 
into account (Kneib et al."1996 |), an d adop ted the lat- 
est lensing model from Richard et al. (2010) to compute 
the lensing magnification. We first derived raw number 
counts directly from the detected sources in our catalogs. 
We then estimated the true number counts and the er- 
rors on the raw number counts using Monte Carlo simula- 
tions. Finally, we calculated the amount of EBL resolved 
based on the final adopted true number counts curves. 
We discuss our procedure in more detail in the following 
subsections. Note that in order to undertake a simple 
and clear analysis, in this section we only used sources 
that were independently and uniformly selected at > 4 cr 
in each bandpass without regard to whether there was a 
submillimeter counterpart in the other bandpass. 

4.1. Delensed Raw Number Counts 

The SCUBA-2 850 /im sources have a coarse resolu- 
tion with a beam size of 14", which can cause large 
uncertainties in estimating the magnifications wh ere the 
sources are strongly amplified ( Chen et al.||201l" ). Thus, 
for 850 fiui sources with a clear counterpart identified in 
either deep radio interferometric maps, 450 /im data, or 
submillimeter interferometric maps from the SMA (see 
Section 5), we adopted the redshifts and the positions of 
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TABLE 1 

SCUBA-2 450 micron detected sources 



Name ID R.A. DEC S/N /450 Counterpart S/Nc,85o /c,850 
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TABLE 5 

SCUBA-2 850 micron detected sources 



Name ID R.A. DEC S/N /gso Counterpart S/Nc,450 /c,450 

(J2000) (J2000) (m.Jy) ID45() (mJy) 
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Fig. 3. — A370 SCUBA-2 S/N maps at 450 /xm (top panel) and 850 /jm (bottom panel) with a linear green scale from -4 to 4. The size 
of both maps is 13' X 13'. Blue squares are the 20 450 um sources given in Table [4] and red circles are the 26 850 fim sources given in 
Table [s] Yellow circles are SCUBA 850 sources from |Cowie et aLl ||2002[ l. White dashed curves are noise contours with levels 0.82 X 
(2, 3) mjy/beam at 850 /xm and 3.92 X (2, 3) mjy/beam at 45U^m. ine sizes of the circles and squares correspond to 2 times the beam 
FWHM of their respective wavelength (Table [Tjl. 
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TABLE 6 

SCUBA-2 850 micron sources compared with the literature 



Source ID 




Chen et al. Ss50 
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TABLE 7 

SCUBA-2 450 micron sources compared with the literature 



Source ID 


S450 


Smail et al. S450 
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Fig. 4. — The result of the completeness tests on A370 at 450 
(asterisks, dashed curve) and 850 /im (triangles, solid curve) vs. 
flux. In the simulation, the flux step at 450 (850) fim is 1 mJy 
(0.5 mJy). For a clear visual demonstration, the spacing between 
the asterisks is 2 mJy, and the spacing between the triangles is 
1 mJy. 



the counterparts for estimating the magnifications. The 
A370 cluster itself is at z = 0.37. For any sources having 
moderate amplifications with redshifts beyond z — 1, the 
amplification is on ly weakly dependent on the redshifts 
(Blain et al.J 1999). Thus, for sources without redshift 
xs, we adopt z = 2.5 based on the current 



measurements, we adopt z 
finding that t he median redshift of the SMG population 
is around 2.5 (Chapman et al. 2005 Wardlow et al. 2011 ). 

Number counts are usually described in two ways: cu- 
mulative number counts or differential number counts. 
When the sample size is small, cumulative number counts 
give a cleaner visual picture of the underlying shape of 
the counts; however, errors on the cumulative counts are 
correlated and need to be carefully taken into account 
when one tries to fit the counts. On the other hand, dif- 
ferential counts have independent errors on each count 
measurement, making it straightforward to do fitting 
on the counts. Our sample size is big enough with at 
least 10 sources in each bandpass that we can focus on 
analyzing the differential counts. In some cases we will 
also present the cumulative counts, but those are based 
on the differential counts. 
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Fig. 5. — At z = 2.5, the source plane areas over which a source 
with a given flux in the source plane would be detected at > 4 cr 
in the image plane. Dashed (solid) curves are for 450 (850) ^m. 
For comparison, we also show the 850 fim image plane area as the 
horizontal dotted segment with the lower flux limit representing 
four times the 850 /im r.m.s. at the center of the image. 



In a blank-field survey, the differential number counts 
can simply be calculated by dividing each source in a 
given flux bin by the area surveyed, summing these mea- 
surements, and then dividing by the flux bin width. How- 
ever, in massive cluster flelds, because of the gravita- 
tional lensing, the area used to divide each source is not 
simply the image plane area. It should be the total source 
plane area in which each source is located and could be 
significantly detected with its flux. 

To flnd the detectable source plane area for each 
source, we gridded the image plane into 1" pixels and 
projected it back to the source plane at the source red- 
shift. The total image plane area is the sum of any part 
of the map that has more than 1/9 of the central ex- 
posure time (Table [2]). We found the spatial boundary 
of the source plane and then calculated the total source 
plane area. The average magniflcation of our observa- 
tions at a given redshift is the ratio between the image 
plane area and the source plane area, which is 1.16 at 
both wavelengths at z = 2.5. We then regridded the 
source plane back into 1" pixels and projected them to 
the image plane. We counted the number of source grid 
points that were detectable at > 4 cr on the image plane, 
assuming each source grid emits the specified flux, and 
then multiplied the total number of points by 1 arcsec^ 
to obtain the total detectable source plane area. In Fig- 
ure [5] we demonstrate the source plane areas over which 
a source with a given flux in the source plane could be 
detected at > 4.0 cr in the image plane, assuming z — 2.5. 

We inverted the source plane area of each individual 
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TABLE 8 

FITS ON THE 450 AND 850 MICRON RAW DIFFERENTIAL NUMBER 
COUNTS 



Wavelengths 


Equation 


No 


So 


a 




450 


(3) 


114 


25.0 


2.10 


7.12 


850 


(3) 


1024 


5.09 


1.51 


4.75 



source to obtain the source surface density. Then, we 
added up the surface density of sources lying within each 
flux bin and divided that sum by the flux bin width to 
obtain the differential counts. The differential counts are 
best described by a broken power law, because once the 
map is sensitive enough to resolve the bulk of the EBL, 
the faint-end slope of the counts needs to be shallower 
than the bright-end slope so that the total EBL is finite. 
We performed fits using the broken power law with 
basically two different slopes on each side of the char- 
acteristic flux. The sources with the characteristic flux 
dominate the contribution to the EBL. The broken power 
law is of the form 




if 5* < 5*0 
if > 5*0 



(3) 



We show the results of our fits in Table |8] We also plot 
them in Figure|6] (black lines), along with the raw number 
counts (black circles). Given the complicated effects of 
gravitational lensing, uneven sensitivity, and confusion, 
the true number counts and the errors on the raw number 
counts are best estimated with Monte Carlo simulations. 

4.2. Simulations 

We first estimated the true noise using the Jackknife 
resampling technique by subtracting two maps that were 
each created from coadding roughly half of the data. 
Since the real sources, regardless of the significance of 
the detection, are subtracted off, the residual maps are 
source-free real noise maps. We then estimated the 
true noise from the noise maps scaled by a factor of 
i/tl X t2/ (tl-|-t2), with tl and t2 representing the ex- 
posure time of each pixel from the two maps. We found 
central noise values of 3.82 and 0.72 mJy/beam at 450 /im 
and 850 /im, respectively, about 0.1 mJy/beam lower 
than the original maps containing faint sources. 

We then took the broken power law fits from Sec- 
tion 4.1 as the underlying models of the source popu- 
lations, drew sources from these models, and populated 
them onto the source plane at z = 2.5 at random posi- 
tions with no clustering. For bright sources, where the 
average number of input sources is less than one, we ran- 
domized their input as either one or zero, with the num- 
ber of each being determined by the probability of seeing 
a source at that flux. We imaged the sources back onto 
the image plane using LENSTOOL, added them into the 
true noise maps, and then ran through the same proce- 
dure of source extraction and number counts calculation 
that we performed on the data. Note that because DIMM 
performs a high-pass filtering in Fourier space to filter out 
the low spatial frequency structures, our data maps are 
not sensitive to the zeroth order scaling factor seen in 
our simulated maps due to many evenly distributed faint 
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Fig. 6. — The raw and simulated differential number counts at 
(a) 450 /im and (b) 850 /tm. In each panel, the black circles are 
the raw number counts, and the black curves are the broken 
power law fits to the raw number counts. The blue dot-dashed 
curves, representing the true number counts, are the final adopted 
model curves for our Monte Carlo simulations (see Section 4.2). 
The red dashed curves and the shaded regions are, respectively, 
the recovered mean number counts and the 90% confidence range 
obtained from the Monte Carlo simulations. The flux ranges en- 
closed by the black dashed lines are where the recovered curves are 
well constrained. 



sources. We therefore applied the filter on the simulated 
maps before we proceeded to extract sources. 

We iterated this procedure 500 times using at least 
10000 simulated data points each time and compared the 
recovered curves from the simulations that connected the 
mean values of each flux bin with the models. We found 
a systematic flux/count boost, which we investigated by 
tracing the input source properties of each detection. For 
each > 4 cr source in the output maps, we identified the 
brightest input source located within the beam area as 
the counterpart and plotted its flux ratio in Figure [Tj 
We did not remove potential false detections if the ra- 
tio were unreasonably high (i.e., a flux ratio > 2), since 
we kept every detected source for the real counts, and 
it is critical to follow the same procedure for the sim- 
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Fig. 7. — The ratio between the delensed fluxes of the detected 
sources and the input fluxes from our Monte Carlo simulations. 
Black dots are ~10000 simulated data points. Red (blue) curves 
are the mean (median) value of the flux ratio in the different S/N 
bins. The area enclosed by the green dashed curves represents the 
1 cr range relative to the mean values. 



ulated counts as for the real counts. At S/N = 4, we 
found median flux boosting factors of 1.61 and 1.31 at 
450 and 850 /im, respectively. At 850 fiui, this is in 
agreem ent with previous SCUBA studies. |Cowie et al.| 
( 2002 1 performed a Monte Carlo simulation with a single 
slope input model of the A370 field and found the boost- 
ing factor to be around 1.25 at 850 ^m. Various SCUBA 
blank-field studies also found similar resu lts ( jEales et aL| 



Scott et al. 112002 Wang et al. 112004 1. 



2000[ 

We then applied the corrections to the models by ad- 
justing the positions of the characteristic fluxes and the 
scaling factors of the broken power laws. Note that we 
kept the measured slopes, since the boost is much more 
sensitive to the S/N than to the intrinsic fluxes. 

We list the parameters of our final adopted model 
curves (hereafter referred to as our true number counts 
curves) in Table [9| and we plot them in Figure [6] (blue 
dot-dashed curves]. We show the recovered mean num- 
ber counts as the red dashed curves. We use gray shad- 



TABLE 9 

The true number counts curves at 450 and 850 micron 
FROM THE Monte Carlo simulations 



Wavelengths 


Equation 


No 


So 


a 


/3 


450 


(3) 


68.4 


20.8 


2.10 


7.12 


850 


(3) 


485 


4.85 


1.51 


4.75 



ing to represent the 90% confidence region. Between 
~ 1.6 — 12 mJy at 850 /xm and ~ 8.5 — 30 mJy at 450 /im 
(black dashed vertical lines), the recovered mean number 
counts curves agree well with the broken power law fits 
to the data (black solid curves) and are well constrained 
with at least 100 data points per bin. The fainter-end 
counts are overestimated by the curves at both wave- 
lengths. Because they are not well sampled due to the 
small source plane area, the output counts are domi- 
nated by serendipitous detections from strongly lensed 
sources, and the number counts curves based on them 
are highly boosted. The same argument can be applied 
to the brighter-end counts, where the curves are more 
unstable and scattered. The surface density of bright 
sources is too low to generate an adequate number of 
sources given the survey area of our observations. Thus, 
the output counts curves are dominated by the boosted 
fainter sources. 

Because our maps have not yet reached the confusion 
limit, the cause of the boost is mostly due to the sta- 
tistical fluctuations of the source flux measurements for 
fl ux- limited obser vations, known as the Eddingto n bias 

f:ddington"1913'). Previous SCUBA studies (e.g., Ealesj 
al.^2000 ) also found that the dominant boosting factor 
is the noise. 

4.3. True Number Counts and the EBL 

We took our > 4 cr sources and deboosted their fluxes 
based on the mean flux boosting ratio shown in Figure [7] 
(red curves). We computed the true number counts in a 
similar way as presented in Section 4.1. However, instead 
of using the delensed fluxes, we added up the surface 
densities of sources with delensed and deboosted fluxes 
lying within each flux bin. Tables [T0| and [TT1 give the true 
differential number counts and the number of sources in 
each flux bin at 450 n m and 850 nm. The errors are based 
on Poisson statistics ( Gehrels||1986 |. The faintest fluxes 
we detected after delensmg and de boosting are ^1.1 mJy 
at 850 /im and ^4.5 mJy at 450 //m. 

We plot our true 450 /im differential number counts 
(black squares) and true number counts curves from our 
Monte Carlo simulations (black curv es) in Figure [8|^a) , 
along with the results from Herschel ( Oliver et al.||2ni0 l 
at 350 /im and 500 um (olive and pink squares, respec- 
tively). In Figure [81 b), we plot our true 850 /im dif- 
ferential number counts (black circles) and true num- 
ber counts curves from our Monte Carlo simulations 
(black curves), along with the results of p revious SCUBA 
850 /im lensing cluster surveys given in Knudsen et al.| 
^08) (purple fille d circles) and o f the SHADES SCUBA 
blank-field survey ( Coppin et al.|2006; blue filled circles). 
We also show the Knudsen et al. (2008) 850 /tm fit to 
their data (blue curve). In Figures Isjc) and (d), we plot 
the cumulative number counts at 450 /tm and 850 /tm, re- 



spectively, which we calculated by summing all the source 
densities above a given fiux. 
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Fig. 8. — SCUBA-2 A370 differential number counts at (a) 450 /im and (b) 850 fim and cumulative number counts at (c) 450 fim and (d) 
850 um . Black filled symbols (squares for 450 /xm and circles for 850 /im) are the delensed and deboosted number countsgiven in Tables [lO] 
and |ll| In all panels, the black curves are our true number counts curves (i.e., the blue dot-dashed curves from FigurelM. Number counts 
obtained from Herschel at 350 fim and 500 /im ( Oliver et al. 2010 1 arc denoted by olive and pink squares, respectively, in panels (a) and 
(c). Blue curves in (b) and (d) show t he 850 Mm fit trom,Knudsen et al.|^2008j . Their number counts are denoted by purple filled circles. 
The blank-field counts from SHADES jCoppin et al.|2006f are denoted by blue filled circles. 



At 450 /Ltm, our bright-end number counts agree with 
the faint-end counts from Herschel (Figure |8]); however, 
our observations probe a factor of ~ 4 deeper. Above 
^ 30 mjy, our curves sUghtly underestimate the Her- 
schel counts, but our survey area is not large enough to 
make this result significant. Because of the steeply rising 
counts as we probe to fainter source fluxes, our 450 fim 
source surface density increases by more than a factor of 
10 compared with the results from Herschel. 

With our 450 fim true number counts curve, the 
amount of 450 fim EBL we directly resolve is ~ 
66.4 Jy/deg2 above 4.5 mJy. This is - 47 - 61% of 
the total 450 /im EEL, depending on the adopted model 

We expect to re- 
mjy based on our 



(Puget et al 
solve 100% o' 



1996 



Fixsen et al.||1998"| 
TEe'i^i} fim EHL at - 1 
counts model, although we note that since our model 



curve does not converge, the curve must turn over to 
a shallower slope in the fainter flux regime. Recently, 
the closest measurements to our results are from the 
deep 500 fim surveys using Herschel. Again, the confu- 
sion limit caused by Herschel' s small aperture size (large 
beam size) hampers its capability to resolve the EBL 
(6%), and the direct counts from Herschel are a factor of 
~4 shallower than the present observations (Figure [S]) . 

A P(D) analysis using pixel flux distributions was used 
to probe deeper in the Herschel maps; however, the re- 
sults were not well cons trained at the fainte r end beyond 



the characteristic flux flGlenn et al. 20101. Bethermin 



et al. (|2012) performed a stacking analysis on a deep 
Herschel 500 /xm map using a sample of 24: fim sources. 
They found they were able to resolve 55% of the EBL 
down to 2 mJy. However, a stacking analysis can only 
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TABLE 10 
450 MICRON True Number Counts 



S450 


N 


dN/dS" 


S45O 




(mJy) 




(mjy-i deg-2) 


(mJy) 


(deg-2) 


5.54 


2 


1455 

-939.7 


4.50 


5786 

-3738 


8.09 


3 


367.5 + 

-200.1 


6.58 


2765 

-1505 


11.83 


2 


83.14 

-53.71 


9.61 


1650 

-1066 


17.29 


3 


161.2 

-87.76 


14.05 


1281 +^^*« 

-697.3 




1 


19.43 

-16.07 


20.53 


236.2 

-195.3 


36.92 


1 


3.76 +^-'' 
-3.11 


30.00 


52.09 +"^-^ 

-43.08 






TABLE 11 






850 MICRON True Number Counts 




N 


dN/dS" 


S850 


N(>S)'' 


(mJy) 




(mjy-^ deg-2) 


(mJy) 


(deg-2) 


1.39 


1 


+14205 

6176 

-5108 


1.10 


7683 +^'«™ 

-6353 


2.13 


2 


553.0 

-357.2 


1.69 


4064 

-2626 


3.27 


8 


1835 +^»^-** 

-634.8 


2.58 


3568 +^^«^ 

-1235 


5.01 


7 


„„. _ +194.8 
361.7 

-133.4 


3.96 


1044 

-384.8 


7.68 


5 


76.44 

-33.02 


6.07 


281.1 +^^"-^ 
-121.4 


11.78 


1 


6.88 

-5.69 


9.30 


34.05 +''*-^^ 

-28.16 



" Differential number counts 
' Cumulative number counts 



provide information on the selected sample as a whole 
and does not reveal the true flux distribution of the indi- 
vidual sources. Our observations directly resolve the con- 
tributing sources, clearly reveal the characteristic flux, 
and provide an unbiased sample of sources for further 
study. 

Moreover, we estimate that without lensing, the de- 
boosted counts can be probed to ~8 mJy at 450 ^m, cor- 
responding to ^30—40% of the 450 EBL. An extremely 
deep SCUBA-2 450 /im blank-field observation can re- 
solve ~70— 90% of the EBL before hitting the confu- 
sion limit (estimated at ~2mJy given the source den- 
sity reaching 1 per 30 beam area) based on our current 
counts model. However, based on the current sensitiv- 
ity of SCUBA-2, taking into account flux boosting, to 
directly detect sources with 2 mJy 450 micron intrinsic 
fluxes at 4 sigma, one needs more than 100 hours of ob- 
serving time, which is not practical. Thus to fully resolve 
450 micron EBL using SCUBA-2, we will need to observe 
the lensing fields. 

At 850 /im, our number counts agree nicely with previ- 
ous SCUBA results in both blank-field and cluster lensed 
surveys. Based on our true number counts curve, our ob- 
servations are able to resolve directly ~ 40 — 57% of the 
850 ^m EBL above 1.1 mJy with the range again corre- 
sponding to the uncertainty in the EBL determination. 
Our 850 /xm observations are not as sensitive as many of 



the confusion-limited blank-field SCUBA surveys; how- 
ever, thanks to the gravitational lensing, we are able to 
resolve more EBL than what has been reported from the 
blank-field SCUBA surveys above 2 mJy (20-30%; e.g., 
Barger et ar]|1999a[ ). 

5. 450-850-RADIO COMPARISON 

The tight, universal correlation between FIR and ra- 
dio lu minosity among no r mal galaxies i n the local Uni- 



verse (Helou et al. 1985 Condon 1992) has been used 



to find the counterparts to the 850 /im SMGs and to es 



timat e their redshifts (e.g., 'Carilli & Yun 1999 Barger 



et al. |2000j). However, this method has not been tested 
on 450 ^m selected sources. With our newly detected 
450 /im sources and the deep 20 cm ma p (1 cr ~ 5.7 /.tJy) 
of the A370 field ( |Wold et ah] 120121 ), we are ideally 
equipped to conduct such a study! Here we study the 
properties of our submillimeter samples in terms of the 
source matches, the positional uncertainties, and the red- 
shift distribution estimated through the various flux ra- 
tios. 

Note that the differential number counts are computed 
by dividing the summed source surface densities within 
a flux bin by the flux bin width, which means that the 
boosting effect on the number counts comes from the 
mean flux boosting factors (red curves in Figure [7]). This 
is shown in Section 4.3 where we used the mean flux 
boosting factors to deboost the raw number counts and 
to compute the true number counts shown in Figure Is] 
that agrees with the underlying true number counts mo<> 
els. 

However, to deboost the fluxes of individual sources 
for characteristic analysis, the median flux ratios (blue 
curves in Figure wh are more representative for the ma- 
jority of the population because of the non-Gaussian dis- 
tribution caused by fake sources with extreme output- 
to-input flux ratios, especially at lower S/N (S/N < 6; 
Figure^. The median flux ratios can be significantly dif- 
ferent than the mean flux ratios, a factor of ~ 2.7 versus 
~ 1.6 at 4 cr at 450 ^m. Thus, for the rest of this paper, 
when we refer to the intrinsic submillimeter fluxes for in- 
dividual sources, they are delensed and deboosted, where 
the latter has been done using the median flux boosting 
factors at the source S/N. 

5.1. Source Matches 

Due to the negative i^-correcti on at 850 ^m and the 
positive JT-correction at 20 cm dCari lli fc Yun| [1999 



Barger et aI1|2000| [Blain et al] [liMj ) , 850 /itm SMUs 



with radio counterparts are biased against high-redshift 
sources, and only a fraction of the 850 /im SMGs are 
seen in the radio. It has been found that around 60% 
of the 850 /xm SMGs with fluxes > 4 mJy have radio 
counterp arts in radio images with a la sensitivity of 



6AtJy (Barger et al.. 
et al.||2007r|Biggs et al." 



2000 



Wang et al. 2004 Ivison 



2011 ) . However, because of the 



I' lR-radio correlation, as the radio images reach deeper, 
there should be a higher fraction of SMGs that have ra- 
dio counterparts. Indeed, the latest study by Barger et 
al. (2012) show that very deep 20 cm interferometric ob- 
servations with a 1 (7 sensitivity of ~ 2.5 fiJy can find all 
SMGs with 850 /im fluxes above 3 mJy 
We cross-identified our su bmillim e ter ca talogs with the 



> 5 (J radio catalog of , Wold et al.| ( |2012[ ), created from 
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a 20 cm survey of A370 with a la sensitivity of ^ 6 /iJy. 
We found that 8 out of 18 (44%) of the 850 sources 
with intrinsic fluxes > 4 mJy have robust radio counter- 
parts within the 850 /xm beam at the > 5 a leveL 

We may reduce our radio S/N threshold if we have 
more accurate submiUimeter positions. We ran a simula- 
tion to estimate the probability of a randomly positioned 
450 /im beam selecting a > 3 ct random noise fluctuation 
on a simulated normal-distributed radio map. We found 
the probability to be ~ 5%. The probability should be 
even lower if there are cou nterparts at other w avelengths 



(e.g., 850 /^m). Note that Biggs et al 



radio sources with 5% of the so-called corrected Poisso 



(2011) considered 



nian probability (Browne & Cohen 1978) to be secure 
counterparts for SiVlGs. Using the 450 /im data, we iden- 
tified four more > 3 cr radio counterparts to the 850 fim 
sources. Thus, in total, 12/18 or 66.7% of the 850 /im 
sources with intrinsic fluxes > 4 mJy have radio matches, 
which is consistent with what has been reported in the 
literature. 

For the 450 /fm sources, the fraction having secure 
(> 3 tr) radio counterparts is 85% (17/20). This high 
match rate could mean that both the 450 /im sources with 
intrinsic fluxes > 6 mJy and the 20 cm emission are pref- 
erentially tracing lower redshift starbursting galaxies, be- 
cause both have positive ivT-corrections at high redshifts. 
Meanwhile, 67% (8/12) of the > 4 cr 450 /im sources are 
detected at 850 /^m, while the recovery rate for > 4 cr 
850 /im sources at 450 /im is slightly lower (13/24; 54%). 

The moderate match between 450 /im and 850 /im sam- 
ple again implies that the SMG population is diverse 
and different selection methods reveal only parts of the 
overall population. For example, SMGs detected only at 
450 /im tend to have a higher dus t temperature and/or 
lower redshift (Casey et al. 2012), whereas the 850 /im 
sources without radio or 450 /im counterparts provide 
the most probable candidate high-redsh ift sources (e.g., 
Riechers et al.||2010| [Walter et al.||2012[ ). 



5.2. Positional Uncertainties 

In the absence of 450 /im observations or submiUime- 
ter or millimeter interferometry, radio sources have been 
widely used to locate counterparts to the S MGs detected 
at 850 /im by s ingle-di sh telescopes (e.g., Barger et al. 



2000 



2005 



[Small et al. 2002; Chapman et al.|2005[|Fope et al. 



This method has been shown to be appropriate 
for many of the sources using high-resolution submiUime- 
ter interferometric observations, though in some cases 
there may be multiple counterparts rather than a sin- 
gle source ([Younger et a l. 2007, 2009] [Wang et al.||2011 



Barger et al.|20l2 " Karim et al.)2012( ) . In contrast to fol- 



lowing up each source individually with interferometers, 
we can statistically examine the method by comparing 
the positional offsets between radio and submiUimeter 
sources from real observations and the offsets obtained 
from Monte Carlo simulations. In theory, if the radio 
emission traces the submiUimeter emission, then the ob- 
servational results should agree statistically with the sim- 
ulations. 

The details of the simulations are described in Sec- 
tion 4.2. We note that in the simulations, at both wave- 
lengths we treated the brightest input source as the coun- 
terpart, and we gridded the maps to 1" x 1" in order to 
match with the real submiUimeter data maps. We plot 
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Fig. 9. — The positional offsets between the output and the 
brightest input sources within the beam area from our Monte Carlo 
simulations (red curves). The areas enclosed by the green dashed 
curves represent the 90% confidence ranges relative to the mean 
values. The diamonds are the offsets between our cataloged (a) 
450 fim sources and (b) 850 /im sources relative to the 5 cr radio 
sources. The triangles in (b) are the offsets between the 450 fim and 
850 /im sources. A370-450.1/A370-850.1, A370-450.18/A370-850.7 
and A370-450.9/A370-450.11/A370-850.5 are plotted with the blue 
diamond in (a), the brown diamond and the magenta triangle in 
(b), respectively. 



our results in Figure[9]at (a) 450 fim and (b) 850 /im. The 
observational offsets oetween the submiUimeter sources 
and their secure radio counterpart candidates are in di- 
amonds in both panels, and the simulated mean offsets 
are denoted by red curves with the 90% confidence range 
enclosed by green dashed curves. The observational re- 
sults generally agree with our simulations, except for two 
sources, A370-450.1/A370-850.1 (blue diamond in Fig- 
urej9];a)) and A370-450.18/A370-850.7 (brown diamond 
in Figure |9][b)). A detailed analysis of these two sources 
is given in the later sections and the Appendix, in which 
we conclude that those two sources are likely to be com- 
plex systems with multiple sources contributing emission 
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at different wavelengths. Taking A370-450.1/A370-850.1 
as an example, recent studies of this source have revealed 
its complex nature, including multiple sources having dif- 
ferent characteristics being located within 1 — 2" of one 
another (more details are given on this source in the Ap- 
pendix) and evidence for a significant offset between the 
dust and radio emission ( Genzel et al.|2003 Ivison et al. 
'2010|). Our statistical diagnosis serves as a simple test 



as to whether there is a common origin for the emission 
seen in the different wavebands. 

We also plot the offsets between the 450 fim and 850 /im 
sources on Figure |9][b) (black triangles). Again, we find 
good statistical agreement between the observational off- 
sets and the offsets that we obtain from the simulations, 
which supports the idea that the 450 /xm sources are a 
good positional tracer for the 850 /im sources. Although 
the 450 /im positions are less accurate than the radio po- 
sitions, the 450 /im emission has the advantage of directly 
tracing the thermal dust emission, just as the 850 /im 
emission does. 

Similar to the situation for A370-450.1 mentioned 
above, there is a complex system in the submillimcter, 
A370-450.9/A370-450.11/A370-850.5, that lies off the 
simulated curve in Figure [9](b) (magenta triangle) . This 
could be another example of a multiple, where differ- 
ent sources contribute emission at different wavelengths, 
and sources sele cted at different w avelengths may have 
different origins (Wang et al.||2011] Barger et al. 2012). 
Note that in this plot wc show the offset between A370- 
450.11 and A370-850.5. The offset between A370-450.9 
and A370-850.5 is even larger. 

In general, there is a slight upward bias for the real 
data relative to the simulated mean curve, which could 
be caused by the fact that some of the counterparts may 
not be correct. As we show in subsequent sections, some 
of the sources have multiple radio counterparts within 
the beam area, which is an indication of a complex na- 
ture. Including incorrect counterparts leads to a large 
dispersion and, in our case, an upward bias. 

5.3. Millimetric Redshifts 

We estimated the redshifts of each source using the 
850 /im-to-450 /im and submillimetcr-to-20 cm flux ratios 
assuming a modified blackbody SED template with /? = 1 
(dust emissivity) and = 47 K (dust t emperature) froni 
the local st arbursting galaxy Arp 220 ( Klaas et al.||1997 
Barger et al. 2000 ) . In Table [121 we show the redshifts 
of each individual source estimated using three differ- 
ent flux ratios. We show dots when there is not a mea- 
surement for both fluxes (e.g., Z4^Q/radio on A370-850.il) 
or when there are multiple counterpart candidates (e.g., 
A370-450. 17/850.4 with its multiple radio counterparts). 
In many cases, the redshift estimates using different flux 
ratios, in particular 850 /im-to-450 /im and 850 /tm-to- 
20 cm, agree with one other to within the uncertainties. 
A noticeable disagreement is seen for SMGs not detected 
at 850 /im (e.g., A370-450.6, 450.8, 450.10, 450.12). Here 
the flux ratios between 850 /im and 450 /im indicate a 
low redshift, while the flux ratios between 450 /im- and 
20 cm indicate higher redshifts. Given that 450 /tm data 
suffer from a positive if-correction, much as 20 cm data 
do, such sources are more likely to be at low redshifts. 
The fact they have abnormally high 450 /im fluxes prob- 
ably indicates that they have higher dust temperatures 




2 3 4 
Redshift 

Fig. 10. — Histograms of the source redshifts from Table [12] The 
redshift distribution for the sources detected at both 450 fim and 
850 /im (green shading) has a median z ~ 2.5. The red histogram 
outhnes the redshift distribution for the 850 /tm sources only, while 
the blue histogram outlines the redshift distribution for the 450 /im 
sources only. The median redshift for all the 450 (850) fim sources 
is z ~ 2.3 (z ~ 3.4). A370-450.13 is not included in this plot, since 
in this case the redshifts estimated from the different flux ratios 
diverge. 



than what we assumed in the SED template. Thus, the 
850 /im fluxes lie below the detection limit. 

In the final column of Table |T2] we give our adopted 
redshift for each source that hadgood constraints on the 
redshift estimation. Redshifts with three digits after the 
decimal point represent the spectroscopic measurements. 
We primarily adopted the redshifts estimated through 
the 850 /im-to-20 cm flux ratio, since it has recently been 
shown to be a reasonable redshift estimator with good 
accuracies (Barger et al. 2012). We adopted the redshifts 
estimated from the 850 /im-to-450 /im flux ratio whenever 
the 850 /im-to-20 cm flux ratio was unavailable. 

We plot histograms of our redshifts in Figure [TO} 
Sources detected at both 450 /im and 850 /im (green 
shading) have a median redshift of z ~ 2.5. The red- 
shift distribution of the sources detected only at 450 
(850) /im is denoted by blue (red) lines. We adopted 
the upper limits of the redshifts estimated using the 
850 /im-to-450 /im flux ratio for the sources detected only 
at 450 /im. For 850 /im sources without good redshift 
constraints, we adopted the lower limits estimated using 
the 850 /im-to-20 cm flux ratios. In the final column of 
Table [T2j we show these redshift limits with parenthesis. 
All the 850 /im sources have a median redshift of z ~ 3.4, 
while all the 450 fim sources have a median redshift of 
z ^ 2.3. These results generally agree with the picture 
that 450 /im sources preferentially trace lower redshift 
dusty sources, while 850 /im sources tend towards higher 
redshifts. We do, however, caution that because there 
are many issues about fitting the FIR SED template for 
SMGs such as the degeneracy betwee n source red shift, 
dust temperature and dust emissivity ( Casey|[2012 ), our 
results are first order estimations and should be used 
carefully. 



6. COUNTERPART IDENTIFICATIONS AND DISCUSSION 
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TABLE 12 
Redshifts 



ID850 



2850/450 



^4S0/radia 



A370-450.1 
A370-450.2 

A370-450.3 

A370-450.4 

A370-450.5 

A370-450.6 

A370-450.7 

A370-450.8 

A370-450.9 

A370-450.10 

A370-450.il 

A370-450.12 

A370-450.13 

A370-450.14 

A370-450.15 

A370-450.16 

A370-450.17 

A370-450.18 

A370-450.19 

A370-450.20 



A37Q-850.1 
A370-850.3 

A370-850.2 

A370-850.18 

A370-850.6 



adopted 



+0.76 
-0.00 
+0.81 
-0.84 
+ 1.73 
-0.00 
+ 1.30 
-1.31 



A370-850.8 2.52 



0.00 
2.08 
0.00 
3.48 ^ 
<0.01 

+ 1.50 



-1.61 



<0.80 



0.97 
1.85 
1.43 
5.95 
5.95 

1.74 



+ 0.07 
-0.07 
+ 0.19 
-0.18 
+ 0.22 
-0.22 
+ 0.00 
-2.30 
+0.00 
-0.98 



+0.36 
-0.33 



0.84 
1.90 

1.18 
4.47 



+0.07 
-0.07 
+0.13 
-0.12 
+0.18 
-0.20 
+ 1.15 



A370-850.5 



A370-850.5 



+ 1.48 
-0.79 



1.90 



4.47 



2.52 



2.800 
1.060 

+0.13 
-0.12 

1.519 

+ 1.15 
-0.76 

(0.01) 

+ 1.50 
-1.61 

(0.80) 
(0.94) 





<1.80 


3.67 


+ 2.28 
-1.13 






(1.80) 


A370-850.26 


0.66 

-Q.66 


>5.95 


>3.22 






A370-850.17 


1.52 

-1.52 


2.99 


+ 1.55 
-0.83 


2.53 

-0.49 


2.53 


+0.58 
-0.49 


A370-850.9 


3.31 

-1.70 


1.95 


+0.48 
-0.43 


2.29 

-0.30 


2.29 


+0.31 
-0.30 


A370-850.10 


3.63 

-1.70 








3.63 


+ 1.70 
-1.70 


A370-850.4 


4.61 

-1.45 








4.61 


+ 1.60 
-1.45 


A370-850.7 


5.01 

-1.52 








5.01 


+ 1.71 
-1.52 


A370-850.25 


2.28 

-2.28 


1.42 


+0.34 
-0.35 


1.58 

-0.27 


1.58 


+0.25 
-0.27 


A370-850.15 


2.80 

-2.38 


3.04 


+ 1.87 
-0.94 


3.08 

-0.55 


3.08 


+0.65 
-0.55 


A370-850.il 


>2.93 






>4.29 




(4.29) 


A370-850.12 


>2.14 






>4.25 




(4.25) 


A370-850.13 


>1.98 






>4.71 




(4.71) 


A37Q-850.14 


>1.99 






2.34 +°'''> 
-0.37 


2.34 


+0.39 
-0.37 


A370-850.16 


>1.78 












A370-850.19 


>0.81 






>4.12 




(4.12) 


A370-850.20 


>0.00 






>4.49 




(4.49) 


A370-850.21 


>1.50 






>4.12 




(4.12) 


A370-850.22 


>1.31 






>3.61 




(3.61) 


A370-850.23 


>1.23 






3.08 

-0.63 


3.08 


+0.79 
-0.63 


A370-850.24 


>0.48 






>3.41 




(3.41) 



"Values with three digits after the decimal point are spectroscopic redshifts. 

In this Section, we restrict our sample to the 14 sources 
that are detected at > 5 cr in at least one submillime- 
ter waveband for their robust detections and better con- 
straints on the astrometry (Section 5.2). We refer to 
these 14 sources as our SCUBA-2 robust sample. We 
note that most part of our 450 fim map has a sensitiv- 
ity lower tha n the the Hers chel confusion limit (la 
6.8 mJy; Nguy en et al.|20l"0 | at similar wavelengths, and 
even at 5 cr level, majority (10 out of 14) of our SCUBA-2 
robust sample would not be found by Herschel at 4cr. 

Except for the Spitzer and Chandra data, the data 
at other wavelengths fully cover the SCUBA-2 surveyed 
area. We show the schematic Spitzer and Chandra cov- 
erage in Figure 11 where the 450 /im and 850 /im sources 
in the SCUBA-2 robust sample are marked with red 
and blue symbols, and the identification numbers of the 
850 /Ltm sources are given. The Spitzer coverage at vari- 
ous wavelengths is denoted by the green lines, while the 
gray area sh ows the area observed by Chandra. 

In Figure 12 we give postage stamp images for each 
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Fig. 11. — The schematic sky coverage of the A370 field. The 
approximate SCUBA-2 coverage is denoted by the red circle with 
a radius of 7'. The Spitzer 3.6 /xm and 4.5 /.im coverage is marked 
with the green solid lines, and the 5.8 /^m and 8.0 /xm coverages are 
marked, respectively, with the green dashed and dotted lines. The 
blue dots (red circles) represent the locations of our 450 (850) /tm 
sources in our SCUBA-2 robust sample. The numbers given are 
the identifications of the 850 fira sources. The gray area shows the 
X-ray coverage obtained by Chandra. 

of the 14 sources, including an inverse gray scale z-band 
image, a false color J, H, image, an inverse gray scale 
4.5 /im image, and an inverse gray scale 20 cm image. We 
use cyan and red circles to represent the positional error 
of the 450 /im and 850 /fm sources with 95% confidence. 

In Table 13 we give the delensed photometry of each 
candidate counterpart from the NIR to the radio. We 
label these candidate counterparts with the submillime- 
ter source identification followed by a small "a" , "b" , or 
"c", if there are multiple candidates. We use dots in the 
table when the submillimeter source lies outside the cov- 
erage of a particular wavelength image. We discuss the 
individual sources in the Appendix. Note that because 
there are no sources in our SCUBA-2 robust sample that 
are only detected at 450 /im, when referring to individual 
sources hereafter, we will use the 450 /im identifications 
for sources detected at both 450 /im and 850 /im, and we 
will use the 850 /im identifications for sources only de- 
tected at 850 /im. 

6.1. X-ray Candidate Counterparts 

While a number of studies have shown that most of 
the FIR luminosity from SMGs comes from star forma- 
tion, the issue of what fraction of SMGs contain an AGN 
is far fro m settled. Results hav e ranged from as many 
as 75% ( jAlexander et al.||2005] radio-detected SMGs) 
down to ---^10-30% (|Georgantopouios et al.|2011 Spitzer- 
detected SMGs). Existing Chandra observations of the 
A370 field partially cover our SCUBA-2 survey area and 
include 9 850 /im sources and 8 450 /im sources from the 
SCUBA-2 robust sample. We c ross-correlated the source 
catalogs in |Barger et al. ( |2001 ) and found 2 counterparts 
in hard X-rays (A370-450T7S370-450.2) and 1 in soft X- 
rays (A370-450.16). This gives a fraction of 3/9 robust 
850 /im sources containing an AGN (^33%), which is 
similar to what Georgantopoulos et al. (2011) found but 
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850.1/450.1 




Fig. 12. — Postage stamp images of the 5cr submillimeter sources in our SCUBA-2 robust sample. The images from left to right in each 
panel are an inverse gray scale 2:-band image, a false color near-Infrared (J, H, K) image, an inverse gray scale 4.5 lira image, and an inverse 
gray scale 20 cm image. The images are plotted with a dynamical range from 1% to 99% using a linear scale. The angular size of each 
image is 22" X 22". Cyan and red circles represent, respectively, the positional uncertainties with 95% confidence at 850 /im and 450 /tm 
(see Section 5.1). 



much lower than the results of Alexander et al. (2005). 
Note that both of these studies from the literature used 
much deeper X-ray maps with exposure times of 1 Ms 
and 2 Ms, respectively. Thus, it is not surprising that we 
find a lower fraction of submillimeter sources containing 
an AGN. For our robust 450 fxm sources, the fraction is 
3/8 (-38%). 

6.2. IR Candidate Counterparts 

The IR properties of the 850 /ini SMGs have been ex- 
tensively st udied from the NIR to SpitzerS.6 /i m to 24 /iir 
( Frayeret a l. 2003, 2004; Clements et al. 2004; jPope et al 



2UU5||Hainiine et al |2009( |Wang et al.|2012p . Une of the 



seek color criteria that could efficiently select the correct 
counterparts to the 850 /im SMGs, especially given their 
large positional uncertainties. 

While the NIR color spans a wide range, and there is 
not yet a single selection method that is successful in 
choosing the SM Gs, a MIR color selection proposed by 
Yun et al. (2008) using the 8.0 /im/4.5 /im versus 5.8 fim/ 



3.fc) /im hux ratio diagram seems to be able to separate 
SMGs from other field galaxies. It has been argued that 
this MIR color selection is too general in that it includes 
more than 50% of the fi eld galaxies in the GOODS-N 
field (Hainline et al. 2009[). However, it is still a good first 
Haa 



goals of studying the IR properties of the SMGs is to 



order test on our candidate counterparts. We plot any 
available candidate counterparts from Table [T3| that have 
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TABLE 13 
NIR, MIR, AND Radio photometry 



Source ID 


R.A 


a 


DEC 


a 


offset'' 


J 


H 


Ks 


3.6 //m 


4.5 /jm 


5.8 /im 


8.0 /im 


20 cm 


z 






(J2000.0) 


(J2000.0) 


(") 


(mag) 


(mag ) 


(mag ) 


(mag) 


(mag) 


(mag) 


(mag) 


(MJy) 






850.1/450.1 


2 39 


52.0 


-1 35 


58.5 


2.0 


21.8 


21.4 


21.2 


20.4 


20.0 


19.4 


18.9 


164.5 


2 


8000 


850.2/450.3 


2 39 


52.1 


-1 32 


07.2 


0.8 


21.8 


21.0 


20.2 


19.2 


18.9 




19.1 


150.4 






850.3/450.2 


2 39 


56.6 


-1 34 


26.6 


1.3 


19.9 


19.7 


19.2 


18.4 


18.4 


18.3 


18.3 


245.5 


1 


0600 


850.4/450.17a 


2 39 


38.3 


-1 31 


17.4 


2.1 


>23.5 


>23.0 


>23.0 


22.8 


22.1 






61.7 






850.4/450.17b 


2 39 


38.0 


-1 31 


16.6 


2.4 


>23.5 


>23.0 


>23.0 


23.4 


22.7 






25.8 






850.5/450.11a 


2 39 


39.0 


-1 38 


25.8 


2.9 


23.2 


22.6 


21.9 




22.0 






275.0 




8219 


850.5/450.11b 


2 39 


38.8 


-1 38 


28.1 


2.0 


23.0 


>23.0 


>23.0 


>23.3 


>23.3 


>23.0 


>22.6 


<16.7 






850.5/450.9 


2 39 


38.5 


-1 38 


14.6 


0.4 


>23.5 


>23.0 


>23.0 




22.4 






18.4 






850.6/450.5 


2 39 


48.8 


-1 37 


53.6 


3.9 


>23.6 


>23.1 


>23.1 


21.9 


21.3 






18.3 






850.7/450.18a 


2 40 


01.8 


-1 35 


34.8 


4.0 


18.9 


18.8 


18.6 


19.2 


19.4 


19.7 


20.2 


<14.4 






850.7/450.18b 


2 40 


01.7 


-1 35 


33.2 


2.2 


>23.7 


>23.2 


>23.2 


>23.5 


>23.2 


>23.2 


>22.8 


23.3 






850.7/450.18C 


2 40 


01.6 


-1 35 


29.7 


2.0 


>23.7 


>23.2 


>23.2 


>23.5 


>23.2 


>23.2 


>22.8 


17.9 






850.8/450.7a 


2 39 


53.1 


-1 37 


47.5 


2.5 


20.2 


20.3 


20.1 


20.7 


20.9 


20.9 


22.1 


<14.8 






850.8/450.7b 


2 39 


52.9 


-1 37 


45.3 


1.3 


>23.7 


>23.2 


>23.2 


22.5* 


21.5* 


22.3* 


21.4* 


23.8 






850.8/450.7C 


2 39 


53.0 


-1 37 


45.9 


0.5 


>23.7 


>23.2 


>23.2 


22.5* 


21.5* 


22.3* 


21.4* 


15.2 






850.9/450.15 


2 39 


46.1 


-1 38 


37.0 


1.1 


23.3 


22.2 


22.2 


21.2 


20.6 






67.6 






850.10/450.16a 


2 39 


56.3 


-1 31 


36.7 


2.5 


16.0 


16.1 


16.3 


16.4 


16.6 




14.3 


938.1 




0286 


850.10/450.16b 


2 39 


56.1 


-1 31 


41.1 


2.9 


>23.7 


>23.2 


>23.2 


>23.5 


>23.2 




>22.8 


50.0 






850.13 


2 39 


39.8 


-1 30 


20.2 


1.7 


>23.5 


>23.0 


>23.0 


22.6 


21.7 




21.2 


23.1 






850.18/450.4 


2 39 


45.6 


-1 32 


43.6 


0.6 


20.9 


20.7 


20.2 


19.5 


19.2 




18.9 


144.3 


1 


5190 



all four Spitzer passband measurements in Figure [T3[ a). 
The enclosed area denoted by the solid lines represents 
the color space of the MIR color selection. We plot the 
normal field galaxies selected from a dee p survey 



of th e GOODS-N field as small gray dots ( Wang et al 
20101. For visual reference, the MIR measurements on 
the SMGs observed with the SMA are plotted as triangles 
for the GOODS-N sample (Barger et al. 2012) and as- 
terisks for t he AzTEC/COSMOS sample ( [Younger etaL 
|2007[ |2009[) . Most of our counterpart candidates, except 
A^-lSOTTa and 450.18a, are located within the SMG 
color regions; thus, A370-450.7a and A370-450.18a could 
be false counterparts based on this selection. 
Another interesting way of selecting SMGs in I R color 



^Positions are extracted from the radio counterpart, if present, and otlierwise from the optical counterpart. For 
sources without any radio or optical counterparts, we used the positions from the images with the best resolution. 
'Offsets are calculated with the 450 /xm position as the reference, except for 850.13. 
'Source blended 

parts or that they are located at lower redshifts (e.g., 
A370-450.16a is at z = 0.0286). While the K,-A.5 
color is an indicator for redshifts (Wang et al. 2012), 
the fact that most of our candidates are located below 
the KIEROs color cut is evidence that sources detected 
at both 450 /im and 850 /im are likely at lower redshifts 
than sources selected purely at 850 /im, which is consis- 
tent with o ur p revious results. 

In Figure 14 ^a), we plot Kg versus the 850/im/20 cm 
flux ratio using sources having unambiguous radio coun- 
terparts in our SCUBA-2 robust sample (filled circles) 
along with the secure GOODS-N SMG sample with SMA 
follow-up from Barger et al. (2012) (triangles). We find a 
correlation between these two parameters using the non- 
parametric Spearman's rank correlation test (p = 0.44, 
15 degree of freedom, P-value= 0.04), where the higher 
the submillimeter to radio flux fraction, the fainter the 
Kg magnitude. We derive a minimum fit (black curve) 
of logio{f850 i^ralho cm) = 0.23 X X,(AB) - 3.24.^ We 
als o plot K-h&nd magnitude versus redshift in Figure 
[T4{^b) using any spectroscopic redshifts available for the 
sources in Figure fl4[a) (including the GOODS-N sample) 
and millimetric redshifts (shown with error bars) for the 
ones without spectroscopic redshifts. While the majority 
of the sources have spectroscopic redshifts, the data dis- 
tribution is very similar in both panels, indicating that 
the 850 /im/20 cm flux ratio is a good indicator of red- 
shift. The correlation we found could mean that the 
bright SMG population has a similar rest-frame optical 
luminosity, and the observed Kg-hmiA dimming is caused 
by the source distance. A similar correlation with K- 
band magnitude versus redshift for SCUBA sources was 
found by Clements et al. (2003) and Pope et al. (2005). 
Our correlation reasonably agrees with their results. 



space 


is to 


use t 


et al. 


20101 


20121 


by the traditional 



-4.5 /im color (KIEROs; Wang 
As opposed to the sources selected 
r ed color criteria, su ch as DRGs (Dis- 



tremely R ed ObjectsT lErston et al. 19881 |Hu fc Ridg^ 
[way' 1994), KIEROs with - 4.5 /im (AB) > 1.6 are 
mostly dusty starbursting galaxies, and one third of the 
SMGs that h ave submillimeter interferometric follow-up 
are KIEROs ( |Wang et al.|[2012| . We examine our coun- 



terpart candidates on a K^-^.b ijwi versus 3.6 /im-4.5 /im 
diagram in Figure |13[ b) . This is similar to Figure 1 of 
Wang et al. (2012). We again plot the normal field galax- 
ies selected from a deep Kg survey of the GOODS-N field 
as small gray dots as in Figure 13 a) . The SMGs in the 
GOODS-N field with robust counterpart identifications 
from Barger et al. (2012) are plotted as triangles, and our 
candidate counterparts are plotted as filled circles. Most 
of the SMGs cluster at the top-right corner, indicating 
a very red nature. Three of our candidate counterparts, 
including A370-450.7a and A370-450.18a, are located be- 
low the SMG locus in ifs-4.5 /im color, which could mean 
that either they are simply not the true SMG counter- 
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Fig. 13. — IR color-color diagram with filled circles representing 
our SCUBA-2 robust sample sources with measurements in all four 
passbands and triangles showing the GOODS-N SMG sample ob- 
served with the SMA by Barger et al. (2012). (a) 8.0 /jm/4.5 /^m 
vs. 5.8 /im/3.6 /im color-color diagram. Another SMG sample ob- 
tained with the AzTEC/COSMOS survey at 890 ^tm and observed 
with the SMA is denoted by asterisks (Younger et al. 2007, 2009). 
The area enclosed by the solid lines is the color space proposed by 
Yun et al. (2008) for selecting SMGs. (b) _ft's-4.5/im vs. 3.6 /xm- 
4.5 /im color-color diagram. The field galaxies selected in a deep 
Ks survey of the GOODS-N field by Wang et al. ( 2012) are plot- 
ted in gray dots, and the area above the horizontal line shows their 
KIEROs selection {Ks - 4.5 /xm (AB) > 1.6). Two sources (450.7a 
and 450.18a) have the sameflux ratios, so we have moved the data 
point of 450.7a 0.05 mag to the left for clarity. 



In this paper, we present results from our SCUBA-2 
450 /im and 850 [ivn observations of the massive lensing 
cluster field A370. Our observations reach a 1 cr sensitiv- 
ity of 3.92 mJy/beam and 0.82 mJy/beam at 450 /im and 
850 [xm, respectively. We summarize our findings below: 

1. We detected in total 20 sources at 450 //m and 26 
sources at 850 /im with a significance of > 4 cr. Using the 
latest lensing model for A370, we derived the intrinsic 
submillimeter fluxes and observed number counts. We 
ran Monte Carlo simulations to obtain the true number 
counts that take into account lensing effects, complete- 
ness, noise distributions, and flux boosting. 

2. We found that our observations resolve ~47-61% 
of the 450 //m EBL above 4.5 mJy and ~40-57% of the 
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Fig. 14. — (a) 850/xm/20 cra-Ks color-magnitude diagram. Our 
sample sources are shown in filled circles with their IDs. The lin- 
ear solid curve represents the minimum x^sqruare fit to our data 
points. Results from submillimeter interferometric follow-up on 
GOODS-N SMGs are plotted in triangles (Barger et al. 2012). (b) 
Redshift— diagram using the same sample sources as described 
in (a). Millimetric redshifts are adopted with error bars whenever 
the spectroscopic redshifts are not available. 



850 /im EBL above 1.1 mJy. The 450 /im counts are the 
deepest determination to date at this wavelength. 

3. We used a statistical method to examine the ori- 
gin of the 450 /im and 850 /im emission by plotting the 
observational positional offsets between the 450 //m and 
850 /zm emission against the ideal distribution of offsets, 
assuming common origins, using our Monte Carlo simu- 
lations. We found statistically that the 450 /xm emission 
is a good positional tracer for the 850 /im population. 

4. We estimated the redshift distributions for our sub- 
millimeter samples using millimetric flux ratios assuming 
a modified blackbody SED template from the local star- 
bursting galaxy Arp 220. We found a median redshift of 
^ 2.5 for only the sources detected at both 450 /im and 
850 /im, z ~ 2.3 for all the 450 /im sources, and z ^ 3.4 
for all the 850 /im sources. 

5. The percentage of SMGs having secure radio coun- 
terparts is 85% for 450 /im sources with intrinsic fluxes 
> 6 mJy and ~67% for 850 /zm sources with intrinsic 
fluxes > 4 mJy. 67% of the > 4 cr 450 /im sources are 
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detected in 850 /ini, while the recovery rate for 850 /im 
sources at 450 /im is ^54%. The recovered rate we found 
is consistent with the scenario that both 450 /xm and 
20 cm emission trace lower redshift dusty starbursting 
galaxies, while 850 /im emission comes from dusty sources 
located at higher redshifts. 

6. We identified potential counterparts in various 
wavelengths from the X-ray to the MIR for our SCUBA- 
2 robust sample (S/N > 5). Three X-ray counterparts 
are found to be correlated with our robust submillimeter 
sample. The AGN fraction of our 450 /zm sample is 3/8, 
38%, and it is 3/9 (33%) for the 850 sample. The IR 
colors of our candidate counterparts mostly agree with 
the locations of SMGs in previous diagrams, except for 
a few sources that may be incorrectly identified. We 
found a correlation between the 850/xm/20 cm flux ratio 
and the /Tg-band. 
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APPENDIX 
INDIVIDUAL SOURCES 

A370-850.1 (450.1) - This is the brightest submiUimeter source at both 450 /im and 850 /im in the A370 region and 
the first submiUimeter source discovered by SCUBA (Smail et aL 1997). Follow-up observations with high spatial 
and spectral resolution have revealed its complex n ature . Based on several optical and CO spectral line observations 
pYayer et al.||1998| |Ivison et al.||1998| [Genzel et al^2003[ ), one of the two optical counterparts, LI, was first identified 
as the main submiUimeter emitter, possibly due to its AGN and star-formi ng activity, while th e nature of the other 
optical source, L2, was debated between scattered AG N light from LI (Vernet & Cimatti 2001) or the remnant of a 
merger (Ivison et al. 1998). The most recent studies by Ivison et al. (|2010j) using the CU J=l— )-0 line and data newly 
obtained from Hubble and Spitzer argue that a heavily dust-obscured starburst between LI and L2, L2SW, is more 
likely to be the true counterpart. They also argue that together with a bright UV emitter, LIN, this could be a merger 
system with at least two galaxies. However, due to the positional uncertainties of all the sources, it is hard to rule 
out the scenario proposed by Genzel et al. (2003) of a gas reservoir residing in a massive, extended disk around LI. 
Interestingly, our 450 /im emission peaks on LI, while the 850 /im emission peaks on L2SW. Based on our results, we 
argue that both LI and L2SW contribute to the FIR emission. Future subarcsecond observations using the EVLA or 
ALMA will be the key to fully understanding a system with such complexity. For simplicity, in this paper we only 
provide the photometric information for LI. 

A370-850.2 (450.3) - With the help of better spatial resolution data at 450 /im, an optically faint, MIR and radio 
bright counterpart is clearly identified. The redshift estimate for this source using the 850 /im-to-20 cm fiux ratio (see 
Section 5.3) is ~ 1.90 ±0.12. 

A370-850.3 (450.2) - This source was also detected in a SCUBA survey at both 450 /im and 850 /im (Smail et al. 
2002). Two possible counterparts are detected across all bandpasses. Detailed observations with the Keck II LRIS 
spectrograph show that L3, an AGN at z = 1.06 with a ring morphology, is more likely to be the counterpart than L5, 



a passive cluster elliptical (Barger et al. 1999b| . This was confirmed by high-resolution interferometric observations 
using the SMA (Chen et aL 2011), and our high S/N 450 /im detection confirms this result. 

A370-850.4 (450.17) - Both a 5cr and a 4cr radio source are located within the submiUimeter positional errors. Both 
sources are optically and IR faint and therefore very dusty. It is likely that both radio sources contribute part of the 
submiUimeter emission. 

A370-850.5 (450.9/450.11) - The submiUimeter morphology of A370-850.5 extends towards A370-450.9. We sub- 
tracted the PSF centered at the peak position of A370-850.5 to obtain the 850 /im flux for A370-450.9, which is 
3.91±1.49 mJy. A 3(T radio source sits close to the A370-450.9 position with faint emission at 4.5 /xm. A bright radio 
source candidate counterpart (A370-450.11a) with a spectroscopic redshift of z = 3.8219 is located on the verge of 
the 95% confidence positional area. Another candidate counterpart (A370-450.11b) is a faint z-band source. This 
region is the most complex system in our survey with potentially up to three sources (A370-450.9, A370-450.il, and 
A370-850.5) contributing to the submiUimeter emission. 

A370-850.6 (450.5) - One Spitzer source is located within the submiUimeter positional error. It has a 3cr radio 
counterpart and is not visible in the optical and NIR. From the 850 /im-to-20 cm fiux ratio, we estimated a redshift of 
z ~ 4.47 for this source. 

A370-850.7 (450.18) - A bright galaxy candidate counterpart (A370-450.18a) is located very close to the 850 /im 
position. This gal axy was reported in earlier optical and NIR cluster survey work ( [Stanford et al.|1995 Margoniner & 
de Carvalho|2000 ). We estimated the redshift of the SMG to be z ^ 5 from the 850 /im-to-450 /im tlux ratio. Although 
there is no redshift measurement for this source, the fact that it is resolved in the optical and IR makes it likely to 
be at low redshifts. In addition, the fact that the SED of A370-450.18a decays towards longer wavelengths makes it 
unlikely to be the origin of the submiUimeter emission. More tentative candidate counterparts are selected within the 
450 /im positional error, where there are two > 3.5 a optically faint, IR faint, radio detections. If those two radio 
sources are proven to be the correct counterparts, then they would be great examples of how 450 /im positions can be 
used to select counterparts to 850 /im sources. 

A370-850.8 (450. 7) - The candidate counterpart A370-450.7a is seen in the optical and IR and was detected in other 
optical {g,r,i) surveys (Margoniner & Carvalho 2000). Faint emission in the Spitzer passbands can be seen towards 
the two faint radio source candidate counterparts, A370-450.7b and A370-450.7c. We estimated a redshift of z ~ 2.5 
for the SMG using the 850 /im-to-450 /im flux ratio. 

A370-850.9 (450.15) - One optically faint but MIR and radio bright counterpart is identified. We estimated a 
redshift of ^ 2.29 ± 0.3 for this source using the 850 /im-to-20 cm fiux ratio. 

A370-850.10 (450.16) - A bright local galaxy candidate counterpart (A370-450.16a) with z = 0.0286 is seen in the 
north -east. This source was reported in var ious observations from the optical to the X-ray (Margoniner & Carvalho 
2000; Barger et al.|200l[ Metcalfe et al.|2003 ) . Although its IR color does not put it around the locus of a typical SMG, 
we cannot rule out the possibility ot A370-450.16a contributing to the submiUimeter emission, since it is bright in the 
radio. Meanwhile, another radio source candidate counterpart (A370-450.16b) is located in the south-west. Its optical 
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and IR faint nature makes it another plausible counterpart. It is likely that both A370-450.16a and A370-450.16b 
contribute part of the submillimeter emission. 

A370-850.18 (450.4) - An optically faint, MIR and radio bright counterpart is clearly identified by the 450 /um 
emission. This source has a redshift of 2; = 1.519 (Wold et al. 2012), which agrees within the errors with our redshift 
estimations for this source using the flux ratios. 



